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^ ; ABSTRACT 

We have modeled the displacement of luminous X-ray binaries from star clusters 
in star-burst galax i es wit h an evolutionary population synthesis code developed by 
■ Hurley et al. ( 200d . 2002 ). In agreement with Kaaret et al. ( 2004 ). we find significant 



o 



o 



spatial offset of X-ray sources from their parent clusters, and the apparent X-ray lu- 
minosity vs. displacement correlation can be roughly reconstructed. The correlation is 
^j"! I not sensitive to the fundamental properties of the clusters (e.g., initial mass functions 

of the binary stars) and the kick velocity imparted to the newborn compact stars, ex- 
cept the common envelope parameter acE- We present the distributions of the main 
■ parameters of the current X-ray binaries, which may be used to constrain the models 

for the formation and evolution of X-ray binaries with future optical observations. 
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'"^ ■ 1 INTRODUCTION 

> J . 11 

, X-ray binaries (XRBs) contribute a significant fraction of the X-ray radiation of normal galaxies ( Fabbiandll989f ). They are 



binary systems containing an accreting neutron star (NS) or black hole (BH) and a normal companion star. Based on the 



00 

, masses A/op of the optical compani ons, XRBs are conventionally div ided into high-mass X-ray binaries (HMXBs) and low- 
mass X-ray binaries (LMXBs) (e.g. IVerbunt fc van den Heuvellll994l ). In HMXBs, the massive (Mop > lOM©) companions 



, generally have strong stellar winds (with mass loss rate ~ 10 ® — 10 ^Mqjx part of which can be captured by the compact 
star; while LMXBs, in which Mop < 1.5Mq, experience mass transfer through Roche-lobe overflow (RLOF) of the companion. 



at a rate of ~ 10 — 10 M© yr . Betw een them are intermed iate-mass X-ray binaries (IMXBs), in which the companion 
stars' masses are in the range ~ 2 — lOM© (|van den Heuvel|[l975l ). Mass transfer in these binaries also occurs through RLOF, 



but on much faster, (sub)thermal timescale of ~ 10'* — 10^ yr. 



' XRBs in galactic disks are thought to have evolved from primordial binaries, in which a high-mass primary star (M > 

^ , IOM0) formed the compact star and the secondary star as its companion. The formation and evolution^ of X RBs are often 



accompanied with mass transfer and loss of mass and orbital angular momentum (ITauris fc van den Heuvej2006l . for a review). 



If the primary star evolves to be a (super)giant and fills its Roche lobe (RL), its mass is transferred to the secondary via 
RLOF. If the mass ratio of the primary and secondary stars is sufficiently high or the primary star has a deep convective 
envelope, the mass transfer process occurs on a dynamical timescale and is highly unstable, so that a common-envelope (CE) 
enshrouding binary results. The secondary star is captured by the expansion of the giant star and is forced to move through 
the giant's envelope. The resulting frictional drag will cause its orbit to shrink rapidly while, at the same time, ejecting the 
envelope before the naked core of the giant star explodes to form the NS/BH. The binary, if it survives the supernova (SN), 
will evolve to be an XRB when the secondary begins to transfer mass to the compact star. Note that during the formation 
and evolution processes of XRBs there may be several instances of mass transfer and CE phases. For example, HMXBs may 
end up in a CE phase, as the NS (or BH) is engulfed by the extended envelope of its companion. If the system survives after 
the CE, an XRB with a Helium companion may be produced. The formation of LMXBs in globular clusters often invokes 
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dynamical process such as tidal capture of a low-mass mai n sequence (MS) star b y a NS IjBailvn fc Grindlavl Il987l ) and 
exchange encounters between an NS and a primordial binary i Davies fc Hansen 19981 ). 

Young star clusters and X-ray sources from them have many interesting aspects on modern astrophysics. Recent data shows 
that compact young massive clus ters contain a rich popu lation of massive stars, evidently following a standard Salpeter-like 
upper initial mass function (IMF: iMassev fc Huntell998l ). However, there are also indicat ions that some compact young mas- 
sive clusters have either a fl atter than normal upper mass function or a cut-off at low mass (|Sternberg||l998l : ISmith fc Gallagher! 
[2001 ; McCr adv etHlboOSl ). Additional ly, observational studies on the mass ratio of binaries have reached widely varying. 



even disparate conclusions (reviewed in (Abtlll983l : iLarsonllioOll ). Recent studies find that there a,re two populations o f sec- 
ondaries in clusters, which leads to a bimodal distribution of mass ratios (for more, see §1.2 in iKobulnickv fc Frverl 120071 . 
and references therei n). For example. Flucv. (,2006.). u sing data from the Ninth Catalogue of Spectroscopic Binary Orbits 
I Pourbaix et al. 2004), reassessed the HogeveenT l 199ol ) study and concluded that the data support an excess of 5 ~ 1 twin 
systems. These cases all reveal that potentially more massive binaries will be produced in these young clusters, which may 
present different observational properties when they turn on X-rays. So investigations on XRBs in young star clusters may 
help explore, besides the formation and evolution of XRBs, the funda mental properti es of the clusters as well as recent star 
formation processes in galaxie s, including stellar population in galaxies (|Liu fc Lill2007l ). massive star formation and evolution 
( Kaper fc van der MeeilboOTi ) 



Using observations from Chandra and NICMOS on board Hubble Space Telescope (Jf^'D. lKaaret et al.l l|2004l ) examined 
the spatial offsets between X-ray point sources and star clusters in three star-burst galaxies. They found that (1) the X-ray 
sources are preferentially located near the star clusters, indicating that the X-ray sources are young objects associated with 
current star formation, and (2) brighter X-ray sources preferentially occur closer to clusters. The displacements of the X-ray 
sources observed in these starburst galaxies are likely due to the motion of the X-ray sources caused by the SN explosions 
and/or dynamical interactions with other stars and binaries in the clusters. The absence of bright X-ray sources with large 
displacements suggests that there is some correlation between the luminosity of an X-ray source and its motion when the X-ray 
luminosity Lx > 10^* ergs~^. They proposed that these high luminosity sources may be RLOF BH-XRBs with intermediate 
mass companions, if emitting isotropically and running away at a speed of u ~ 10 kms~^. If u ~ 50 kms^^ instead, the short 
lifetime (~ 4 Myr) of these luminous sources needs very massive companion stars or alternatively the luminosity of the sources 
decreases with age, which can be examined by an evolutionary population synthesis (EPS) calculation. They also pointed 



out that the correlation appears inconsistent with the highly beamed X-ray emission model (|King et al.ll200ll : iKording et al 



2OO2I : iKaaret et al.ll2003l ) because the delay time between the formation of a BH and the onset of the thermal time-scale mass 
transfer is long enough for the X-ray source to move to ~ 1 kpc from its point of origin. 

The spatial offset between the X-ray point sources and the star clusters, especially the X-ray luminosity versus displace- 
ment correlation is determined by the velocity of the binary after the birth of the NS/BH, the time passed since the SN, and 
the mass transfer process. Spatial dist ribution of X-ray sources in g alact ic environment has boon investigated both observa- 
tionally and theoretically. For example. IVan Paradiis fc White! (119951 ) and lWhite fc Van Parad iis (1996) have investigated the 
spatial distribution of NS and B H LMXBs i n our Galaxy, from which they suggested that the compact objects had received 
a kick during t he SN explosion . Paczvhski ( 1990l ) also studied the spatial distribution of Galactic NSs in the scope of SN 
kicks. Recently Zuo et all 1I2OO8I) mod e led th e spatial distribution of Galactic XRBs, incorporating the kinematic evolution of 



kicked binary systems. Kiel fc Hurley ( 20091 ) calculated the scale-heights as well as the radial and space velocity distributions 
of pulsars, considering their kinematic evolution within the Galactic potential. Howeve r, for XRBs born in star clusters. 



besides mass ejection from the binary (iNelenians et al 



, , I999I: Ivan den Heuvel et al.ll2000l ) and the kick velocity imparted on 

the newborn NS caused by the SN explosion jLvne fc Lorimer 19941). an ejection speed via dynamical interactions in clusters 
i Phinnev fc Sigurdssonlll99ll : iKulkarni et al.|[l993l : ISigurdsson fc Hernauistlll993l ) should also be considered. In our calculation 



we on ly consider the former two mechanisms on the motion of XRBs (i.e., primordial binaries; IWebbink et~aLlll983l : IWebbinkI 
1992[ ) though the third one can also influence the binary population, but not significantly in our situation, as explained below 



Dynamical interactions include tidal capture, physical coll isions and exchange en counters. The formation rate of NS-XRBs 
produced through tidal capture process using Eq. (3) from lVerbunt fc Hud (|1987l ) can be estimated as: 



R = nnanvrcify — 6 X 10 



Ml + M2 d lOkms" 



102pc-3 IQ-ipc- 



Mp 



Rr. 



Vrcl 



-yr pc 



(1) 



where rins and n are the number densities of NSs and other stars, respective ly, Vj-d is the relative velocity between the stars 



at infinity, a ~ ■Kd[2G{Mi + M2)/Vj.gi\ (i.e., Eq. (2) in lVerbunt fc Hutlll987l ) is the cross section of the two passengers with 



small relative velocities, Mi the NS mass, M2 the mass of the companion star, and d is the distance of closest approach 
of the two passengers . Here we adopt (Mi -I- M2) ~ 2M© and dmax ~ lOi?© which is estimated roughly from Eq. (1) in 
Verbunt fc HutI ( 19871 ). Given the typical values of youn g massive clusters as Urei ~ 10 kms~^, the core radius rc ~ 1 pc, the 



central density po ~ W'^Mp , pc"^ (jMcCradv et aLlbooj ). we can get n ^ 2 x 10 



op 

IMF of iKroupa et all (|l993l ). following the approach (i.e., §3.1) of I Verbunt fc HutI l|l987l ). here / ~ 0.2 is the fraction of NSs 



pc_ 



remaining in the cluster. Then we can estimate the predicted number of XRBs in one cluster is A^'x — ^Trr^R x T ~ 10 with 
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an X-ray lifetime T = IQ Myr. So the expected number of XRBs produced through this chann el in the 



Kaaret et al 



sample is very small 1). Exchange encounters can also occur on time-scale of a few Myr (jPortegies Zwart et al. 



mm 



altho ugh the formation rate of XRBs through exchange collisions barely competes with two-body encounters IjHut fc Verbunt 
19831 ). However, during the exchange encounters a lower mass binary star tend to be replaced by a more massive participant. 



Hence the effect of exchange encounters is to modify the mass-ratio distribution of the binaries, making equal-mass binaries 
more likely. We examined this effect in our models (i.e., models M3 and M7) and found no significant changes. 

In the present work, we investigated the kinematic consequences of XRBs from star clusters in star-burst galaxies from a 
theoretical point of view. We used an EPS code to calculat e the e xpected cumulative distribution of XRB displacements from 
their parent clusters. Following the approach of Zuo et al. ( 20081 ). we calculated the spatial offset distribution of XRBs with 
luminosities > 10^® ergs~^. We mainly examined several parameters, such as IMF, the secondary star mass function, common 
envelop efficiency and kick velocity, which may affect the formation, evolution and motion of massive XRBs significantly. 
The objective of this study is to use the apparent X-ray luminosity versus displacement correlation to constrain the model 
parameters and the fundamental properties of clusters in star-burst galaxies. We also aim to explore why such correlation 
exists, which may help understand the nature of the sources and may be testified by future observations. 

This paper is organized as follows. In §2 we describe the population synthesis method and the input physics for XRBs in 
our model. The calculated results are presented in §3. Our discussion and conclusions are in §4. 



2 MODEL 



2.1 Assumptions and input parameters 

2.1.1 binary evolution 

To follow the evolution of XRBs, we have used the EPS code developed bv lHurlev et al. I (|2000l . l2002h in our calculations. This 
code incorporates evolution of single stars with binary-star interactions, such as mass transfer, mass accretion, CE evolution, 
SN kicks, tidal friction and angular momentu m loss mechanics (i.e., mass loss, magnetic braking and gravitational radiation). 
We also updated the original code according to lLiu fc Lil (|2007l l and lZuo et al.l (|2008 ). The values of other adopted parameters 
are the sa me as the default pa rameters in Hurley et al. (2002) if not mentioned. 

In the lKaaret et al.l (|2004l l samples, the ages of the star clusters range from 1 to 20 Myr, we then adopt a constant star 
formation rate for 20 Myfl For each model, we evolve lO'' primordial system^, all of which are init ially binary systems . We 
set up the same grid of initial parameters (primary mass, secondary mass and orbital separation) as lHurlev et al.l ( 20021 ) did 
and evolve each binary on the grid. In the following we describe the assumptions and input parameters in our control model 
(i.e.. Model Ml, listed in Table 1). 

(1) initial parameters 

We adopt the IMF of IXroupa et all l|l993h for the distribution of the primary mass (Mi). For the secondary's mass (M2), 
we assume a uniform distribution of the mass ratio M2/M1 between and 1. A uniform distribution is also taken for the 
logarithm of the orbital separation In a. 

(2) CE evolution 

When mass transfer becomes dynamically u nstable, a CE w i ll engulf the binary . An important parameter in the evolution 
of close binaries is the CE parameter acE ( Paczvhski 19761 : lb en fc Livio 19931 ). It describes the efficiency of ejecting the 
envelope (Afcnv) by converting orbital energy (i?orb) into the kinetic energy that provides the outward motion o f the envelope 
(Eeriy), described as 
199^ 



Webbink 



1984 



QCEAi?orb. There are several schemes for the CE evolution in the literature (e.g.. llben fc Livio 



bv lKiel fc HurlevI l|2006l ) 



Hurley et al.ll2002l ). Here we adopt the estimation of the reduction of the orbital separation suggested 



■ acE 



which yields the ratio of final (post-CE) to initial (pre-CE) orbital separations as 

Of _ McM2 1 
^ ^ Ml McMa/Mi + 2Menv/(QCEArLi) ' 

In the above equations the subscripts / and i denote the final and initial values, respectively; Mc = Afi — Afcnv is the core 
mass of the primary star (Mi) that fills out its RL, tli ~ Rhi/ii is the dimensionless RL radius, and A is a parameter which 
depends on the stellar mass-density distribution, respectively. Here we adopt qce = 1 and A = 0.5. The orbital separation of 
the surviving binaries is often reduced by a factor of ~ 100 as a result of the spiral-in. If there is not enough orbital energy 



(2) 



(3) 



^ We also extended the star formation history (SFH) to 50 Myr, and found very small dlfEerence in the final results. 

^ We also varied the number of the binary systems by a factor of two, and found no significant difference in the final results. 
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available to eject the envelope, the orbital motion of the companion during the spiral-in process may be unable to drive off 
the envelope of the primary star, resulting in coalescence rather a compact binary. 
(3) SN kick velocity 



When a binary survives a SN explosion, it receives a velocity kick due to any asymmetry in the explosion (ILvne &i Lorimei 
1994h . The kick velocity Hk is assumed to be imparted on the newborn NS with the Maxwellian distribution 



exp(- 



Kick' 



(4) 



and we adopt crkick= 190 kms"^ iHansen et al.ll 19971 ') in our control model. For BH systems, we ass ume that t hose BHs formed 
with successful SN explosions (the Carbon/oxygen core mass less than ~ 7.6Mq, Fryer 1999; Fr ver fc Ka logcra were 
imparted on a natal kick velocity, which is inversely proportional to the BH mass, i.e., iikick.BH = 1-4/A/bh x iikick.NS- Otherwise 
no kick velocity is adopted. The velocity (vg ) of the bin ary system is related to both the kick velocity and the orbital velocity 
of the system, and can be expressed as fsee lHurlev et al.ii200^ . for details). 



Ml 



-Vk 



AA/1M2 



(5) 



where Mb = Mi + M2 and — Mb — AMi are the total masses of the system before and after the SN, respectively; 
Mj = Ml — AMi is the current mass of t he primary s t ar aft er losing mass AMi during the SN, v is the relative orbital 
velocity of the stars (expressed as Eq. Al in Hurlev et al. 1 20021 ')). Tidal effect is taken into account to remove any eccentricity 
induced in a post-SN binary prior to the onset of mass transfer. 

We also construct several other models by changing the key input parameters (list ed in Table 1). Recent observations 
indicate that some compact young m assive clusters c ontai n relatively more massive stars ( Sternber3l 19981 : Smith fc Gallagher 
200 ll ). we then consider an IMF from lBallero et al.l l|2007l . i e.. Model M2 and M6), which is more skewed towards high mass 
than in the solar neighbourhoo d for comparison. The secondar y star mass function can also strongly affect the formation 
rate of different type of sources jKalogera fc Webbink|[l998.. : ,Belczvnskil booj ) . For example, diffe rent choice of the secondary 
star mass function can change the formation rates of NS binaries by nearly a factor of ~ 100 ( Frver et al. 19991 ). Here we 
assume the secondary mass following a power-law distribution: P{q) oc q", where q = A/2/M1. Most population synthesis 
studies adopt a flat mass spectrum (i.e., a = in our control model Ml) for systems that are likely to interact, while recent 
observati ons seem to be more consistent with "twins" being a general feature of the close-binary population (i.e., Model M3 
and M7. iDalton fc Sarazinlll995l : iKobulnickv fc Frven 120071 ). As stated before, binary interactions may also tend to modify 
the mass ratio to approach unity. So we also adopt a = 1 to examine its effect. 

■Variations in the CE parameter can considerably change the relative numbers of XRBs. However, reliable value for acE 
has proven difficult to estimate, due to lack of understanding the complicated processes involved. Generally acE — 1 is used 
but its range can change widely from ~ 0.1 to ~ 3.0 (see discussion in §4). So we also adopt acE = 3.0 for comparison (i.e., 
models M5, M6, M7 and M8). 

Finally, since the kick velocity Vk can affect not only the global veloc ity of the binary s ystem but the outcome of the 
XRB evolution, we also adopt akick = 270 kms"^ (i.e., models M4 and M8. iHobbs et al.ll2005i ) for comparison. 



2.1.2 binary motion 

Since star clusters in star-burst regions are usually centrally concentrated, we assume a spherical potential and adopt the 
cylindrical coordinate system (r, 0, z) centered at the cluster's center. The potential of a cluster can be constructed as 

-GM 



$N(r, 2) = 



+ + 



(6) 



where G is the gravitational constant, h the half light radius and M the total mass of stars within the half light radius. Here 

"16 , ,-, r-. 



we adopt M — 1.0 x 10** A/©, and ft = 3 pc ( Ho fc Filippenkol 1996al lbl) in our calculation^. We assume that stars are born 
uniformly in the star cluster with random direction of the initial velocity vector, which gives the initial velocity vectors Vr, 
v^, Wz. Due to the cylindrical symmetry of the potential, two space coor dinates r and z a re sufficient to describe the XRB 
distributions. Then we integrate the motion equations (i.e., Eqs. (19a,b) in |Paczvhskiiri99Gl ) with a fourth-order Runge-Kutta 
method to calculate the trajectories of the binary systems and collect the parameters of current XRBs if turning on X-rays. 
Finally we project the positions of XRBs on the tf) — plane to get the projected distances of XRBs from star clusters, i.e., 
R — ((rcosiys)^ -I- z'^y^'^ where 93 is uniformly distributed between and 2-k. In our calculations, the accuracy of integral is 
set to be 10"'' and controlled by the energy integral. 



We also adopted M = 5.0 X 10^ Mq, and found very small difference in the final results. 
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2.2 X-ray luminosity and source type 



Compact stars in XRBs are powered by eitlier disk or wind accretion. When a star expands to fill its RL as a result of stellar 
evolution, or by angular momentum losses causing contraction of the orbit, it can transfer masses via an accretion disk which 
is fed by RLOF; otherwise, mass-transfer cannot occur except that there is a stellar wind to power an observable X-ray source 
(i.e., wind - accret ion). We adopt the same procedures to calculate X-ray luminosity and divide different types of sources as in 



Zuo et al. 1 2008^ if not mentioned. In particular, for disk-fed sources, the simulated X-ray luminosity form is as follows 



J^boi^outiEdd transients in outbursts 

J7boiniin(Lboi,77EddiEdd) persistent systems. 



(7) 



where ?7boi = 0.3 is the bolometric correction facto r which is introduced to c onvert the bolometric luminosity (Lboi) to the 
X-ray luminosity in the 2 — 10 keV energy range jBelczvnski fc TaamI [iooi l : Lhoi ~ V^a-ccC^ where rj is the efSciency for 



energy conversion, Mace is the average mass accretion rate f or accreting systems and c is the velocity of light; ?7Edd is the 
"Begelman factor" ( Rappaport. Podsiadlowski fc PfahJ 12004 ) to allow super-Eddington luminosities; the critical Eddington 
luminosity LEdd — 'inGMimpc/aT = 1.3 x 10''*(Mi/Mq) ergs~^ (where ar is the Thomson cross section, and rup the proton 
mass). For transient sources, the luminosities in outbursts are taken to be a fraction {r]out) of the Eddington luminosity. Here 
we adopt r^out = 0.1 and 1 for NS systems if the orbital period Porb is less and longer th an 1 day, respectively; for BH systems, 
we adopt 77out — Porb/24 hr and let the maximum peak luminosity not exceed 5I/Edd ( Belczvnski fc TaamI 2003|; Ic hen et al 



199?! : ICarcia et al.|[20oi ). To discriminate transient and persistent sources, we use the criteria of IVan Paradiis (199Q1 for MS 
and red-giant donors, and of Ivanova fc Kalogera ( 20061 ) for white dwarf donors, respectively. 



3 RESULTS 



The three star-burst galaxies M82, NGC 1569 and NGC 5253 in lKaaret et al.l (|200J) all contain several young, luminous star 
clusters and luminous X-ray sources in them. The number of X-ray sources in M82/NGC 1569/NGC 5253 is 42/14/10 with 
the corresponding cluster number 50/58/13, respectively. The spatial offset between the X-ray point sources and the star 
clusters shows that, while the X-ray sources are generally located near the star clusters, brighter X-ray sources preferentially 
occur closer to the clusters, and there is an absence of very bright sources (Lx > 10^* ergs~^) at relatively large displacements 
(> 200 pc) from the clusters. Here we modeled the kinematic evolution of XRBs from the star clusters. The calculated results 
are presented below. 

As stated before we constructed 8 models to investigate how the final results are influenced by the adopted parameters. 
Speciflcally the input parameters in our control model (i.e., model Ml) are SFH= 20 Myr, (Jkick ~ 190 kms"\ a = 0, 
acE ~ 1.0, and the Kroupa IMF. In other models we change only one parameter each time, and the model parameters are 
listed in Table 1. Figures 1 and 2 show the simulated cumulative distribution of the X-ray source displacements (top: Total; 
middle: NS-XRBs; bottom: BH-XRBs) at the age of 20 Myr for models M1-M4 and M5-M8, r espectively. Note th at we only 



select sources in the luminosity range lO'^'* < Lx < 10^* ergs in order to compare with iKaaret et al.l (1200J). We also 



normalize each histogram by t he total number of X RBs within 1 kpc of a star cluster for each galaxy. We find that the results 



are in general agreement with lKaaret et al.l (|2004l . Fig. 2 ). 

Figures 3 and 4 show the modeled distributions of the X-ray luminosities (I/x) at different displacement (i?) from the 
star cluster for models M1-M4 and M5-M8, respectively. The top, middle, and bottom panels are for all XRBs, NS-XRBs, 
and BH-XRBs, respectively. The color bar represents the normalized number ratio of XRBs in the R — Lx plane. Note that 
NS-XRBs have relatively lower maximum luminosities than BH-XRBs in all of the models, due to a lower Eddington accretion 
rate limit. The predicted Lx vs. R relations in all the models are roughly compatible with the observations, but differences 
also exist. For models M1-M4, there is a scarcity of high-luminosity {Lx > 10"^* ergs~^) sources in the 30 — 100 pc region, 
while there is an overabundance of very luminous sources (Lx > 10"^^ ergs~^) in the same region for models M5-M8. In 
addition, the correlations originate from different stellar populations in different models. On one hand, models M1-M4 have 
similar Lx vs. 7? correlations which are constructed by both BH-XRBs and NS-XRBs. Note that BH-XRBs dominate at the 
high-luminosity (Lx > 10^* ergs""*^), small-offset (10 < R < 300 pc) region (named as region A), while NS-XRBs dominate 
at the low- luminosity (Lx < 10^* ergs~^), small-offset (10 < J? < 300 pc) region (named as regions B) and large-offset 
(300 < R < 1000 pc) region (named as region C). On the other hand, for models M5-M8, the correlations mainly result from 
BH-XRBs instead. It is interesting to note that there is a fiat-roofed edge in the 10 — 100 pc region, which is not caused by 
the Eddington luminosity limit for BH-XRBs. 

In order to explore the nature of the XRBs in different regions, we need to examine their observational properties (i.e., 
current mass M2 and spectral type of the donor star, orbital period Porb, and system velocity distribution). We present the 
Lx — M2 (left), Lx — Porb (middle), and Porb — M2 (right) distributions of XRBs in the 10 < P < 300 pc region for models 
M1-M4 and M5-M8 in Figures 5 and 6, respectively. Figures 7 and 8 show the distributions of the same parameters in the 
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region of 300 < J? < 1000 pc. The detailed source types in regions A, B and C for different models are also listed in Tables 2-4, 
respectively. 

Figures 5 and 6 show that, the XRBs in region A are in short orbital periods (~ 1 hr), with relatively low-mass (~ 1-3 
Mq) donors for all models. They are mainly RLOF BH-XRBs with Helium main-sequence (HeMS) companions (see Table 2). 
The binary velocities are ~ 10 — 20 kms^^ for models M1-M4, and ~ 30 — 100 kms^^ for models M5-M8. Considering that 
they have similar evolutionary timescales (see Fig. 9 below), the difference in the velocities explains the different maximum 
offsets of these high-luminosity sources in different models. The XRBs in region B have larger orbital periods and companion 
masses. For example, the orbital periods are ~ 1 — 20 hr for models Ml and M3-M8, and can reach ~ 10^ hr for model M2; the 
companion masses are ~ 1 — 4Mq for models Ml and M4, ~ 1 — 2OM0 for models M5-M8, and can reach ~ 6OM0 for models 
M2 and M3. Table 3 indicates that the majority (> 80%) of these sources are NS-XRBs. They mainly have HeMS companions 
(~ 50% in model M4 to ~ 85% in model M7), with mass transfer either through RLOF or by wind capture. Figures 7 and 
8 show that, the XRBs in region C all have relatively low-mass (~ 1-4 Mq) companions. Their orbital periods are around 
several hours for all models. The typical velocities are ~ 150 — 300 kms~^, much larger than those of high-luminosity sources. 

We tentatively conclude that the difference in the luminosities for the sources in regions A, B and C is caused by different 
types of the compact objects (BH-XRBs vs. NS-XRBs), of the donors (HeMS star vs. MS star), and of the accretion modes 
(RLOF vs. wind-capture). Sources in region A are mainly BH-XRBs with RLOF mass transfer from an HeMS companion, 
hence reaching very high luminosities. Sources in region B have similar properties as those in region A, but are mainly NS- 
XRBs, with relatively lower luminosities. Sources in region C are either NS-XRBs (in models M1-M4), or contain a portion 
of BH-XRBs with MS companions, so they also cannot reach luminosities as high as those in region A. 

The displacement of the XRBs from the star clusters depends on their velocities at the moment of the SN explosion and 
the delay time from the SN to the beginning of RLOF. In Fig. 9 we present the distribution of the delay times for sources in 
regions A (dash-dot-dotted line), B (solid line) and C (dotted line), respectively. We normalize each the histograms by the 
total number of X-ray sources in each region. They generally range from 1 to 8 Myr, but the times are a bit shorter in models 
M1-M4 than in M5-M8. In the latter they peak at ~ 5 Myr. 

Note that the Lx vs. R relation has different origin between models M1-M4 (acE = 1.0) and models M5-M8 (acB = 3.0). 
In the acE = 1-0 cases, it is constructed by both high-luminosity (Lx > 10^* crgs"^) BH-XRBs and low-luminosity (Lx < 10'''** 
ergs~^) NS-XRBs, while in the qce = 3.0 cases, it results from BH-XRBs instead, so we need to discuss them separately. In 
models M1-M4, luminous BH-XRBs are constrained within ~ 300 pc because of their low velocities, while NS-XRBs, which 
arc relatively dimmer and faster, can move to region C within ^ 10 Myr (sec Table 4). In models M5-M8 , a number of RLOF 
BH-XRBs with MS companions are produced, significantly different from models M1-M4. These XRBs, occupying region C, 
have relatively lower luminosities and higher velocities compared to the BH-XRBs with HeMS donors in region A. 

In order to understand the difference in the binary velocities in different models, we examine the distribution of the 
companion masses M2,sNe, and orbital periods Porb.SNe at the moment of SN explosions. It is clearly seen from Eq. (5) that 
the velocity of an XRB is determined by the kick velocity, total mass and ejected mass, orbital velocity (or orbital period) 
when SN occurs. Figures 10 and 11 show the Porb.sNc — M2,snc distribution in regions A. B and C for models M1-M4 and 
M5-M8, respectively. High-luminosity sources in region A (top panel) generally have longer-period (hence smaller orbital 
velocity) and more massive companions, resulting in smaller system velocity (hence smaller offset from the parent cluster) 
than low-luminosity sources (bottom panel). This fact implies that the orbital evolution plays an important role in kinematic 
motion and the spatial distribution of XRBs. 

The formation of XRBs usually invokes at least one CE phase. So the CE parameter acE can influence the Lx vs. R 
relation by affecting the binary orbit distribution after the CE evolution. This further determines the evolutionary state of the 
donor star during the XRB phase and the system velocity after the SN explosion. It has two contrary effects on the formation 
and evolution of XRBs. Larger values of acE can prevent coalescence of a BH/NS and the companion in a compact binary 
during the unstable mass transfer processes, in favor of the formation of compact XRBs; while for some initially wide binary 
systems, a higher value of acE can cause the RLOF mass transfer to a NS/BH not to occur within 20 Myr, decreasing the 
formation rate of XRBs. So variation of acE may lead to different types of XRB populations. In addition, it can affect the 
orbital velocity, and the velocity of the binary. Larger acE means that more energy was used to drive off the envelop, leading 
to wider orbital separation, and hence smaller orbital velocity after the SN explosion. 

To explore the influence of acE in detail, we present an example evolutionary sequence for Mi , M2 , Porb , Lx of an XRB 
source in region C in Figure 12. We consider a primordial binary system in a 448.332 i?© orbit. The initial stellar masses 
are 14.309 and 3.833 Mq for the primary and secondary, respectively. In model M5, the primary first evolves across the 
Hertzsprung Gap (HG), expands and fills its RL at the time 13.7057 Myr, then transfers mass to the secondary star, which 
is still on MS. Due to the large mass ratio, the secondary is unable to accept the overflowing material, and the system enters 
the CE stage. The secondary spirals within the envelope and drives it off on a dynamical timescale, leaving a compact binary 
(a = 14.845Ji0) with the He core (of mass of 3.482M0) of the primary. At the time of 15.7485 Myr, the He star evolves across 
the HG and fills its RL. The system enters the second CE stage, which makes the orbital separation to shrink again. In a 
very short time, the SN explosion occurs and a NS (of mass 1.362M0) is born. The post-SN binary gets a global velocity 
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of ^ 200 kms ^ due to the large orbital velocity before the SN. Now the binary separation becomes 5.395-R0. Because of 
the small separation, the secondary star will fill its RL immediately and transfer mass to the NS, leading to the formation 
of an NS-XRB. At the time of 16.0062 Myr, the NS collapses into a BH with mass of S.OOSMq, and a BH-XRB with an 
MS companion is formed after twice CE evolutions. The X-ray luminosity Lx first rises to ~ 10^^ ergs~^ and then decreases 
gradually with time. At the time of 17.69 Myr, Lx decreases to be ~ 10^* ergs~^ and the orbital period evolves to be ~ 10 hr. 
By now the binary has moved for about 2 Myr since the SN explosion, with a displacement of about 400 pc. The above story 
explains how BH-XRBs of this type (RLOF and MS companion) can move to region C. However, when we decrease the acE 
to be 1.0, there isn't sufficient energy to drive off the entire envelope of the primary star during the first CE phase, leading 
to coalescence of the giant core and the MS star. So no such type of XRBs will be produced. 

Note that in the above case, most of the low-luminosity {Lx < 10'^* ergs~^) BH-XRBs in region C are formed from 
accretion-induced collapse (AIC) of NS systems. Whether AIC of NSs really happens is still under debate, since quite a few 
binary millisecond pulsars , though thought to have experienced extensive mass accretion, seem to have masses not far from 
1.4Mq ( |Bassa et al. |2006| , and references therein), suggesting significant mass loss during accretion. So we examine this effect 
in Figs. 13 and 14, which are the same as Figs. 3 and 4, except that we do not consider the NS— >BH AIC formation channel. 
It is clearly seen that, in the case of acE = 3.0 (i.e., models M5-M8), the Lx vs. R correlation barely exists, but is constructed 
by NS-XRBs instead, and BH-XRBs disappear in region C. We propose that precise velocity measurement of BH-XRBs may 
serve as a possible way to discriminate different BH formation channels. 

To explore the evolution of the Lx vs. R relation, we plot Lx vs. R in Fig. 15 at the time 20, 30, and 40 Myr for models 
Ml (top) and M5 (bottom). Note that the correlation in both cases exists at 20 Myr after the star formation, and disappears 
gradually as time goes on. If we shorten the SFH to 5 Myr, the correlation disappears much earlier, at the time of ~ 25 Myr. 
So the Lx vs. R correlation cannot hold all the time for individual cluster. 



4 DISCUSSION AND SUMMARY 



We have used an EPS code to calculate the spatial offset distribution of XRBs from their parent star cluster in star-burst 
galaxies. We used the apparent X-ray luminosity versus displacement correlation to constrain models of XRBs. Our study 
shows that the correlation can be roughly reproduced with all models considered, but significant differences exist when 
changing the common envelope parameter acE- In the qce = 1.0 cases (models M1-M4), the Lx vs. R relation is constructed 
by both high-luminosity {Lx > 10^* ergs"^), small-offset (10 < i? < 100 pc) BH-XRBs and low-luminosity {Lx < 10^* 
ergs"^), large-offset (300 < R < 1000 pc) NS-XRBs, while for qce = 3.0 (models M5-M8), it is mainly determined by 
BH-XRBs. In this case, a bunch of high-speed (~ 150 kms~^), low-luminosity, RLOF XRBs with MS donors dominate the 
sources in the 300 < R < 1000 pc region, which are not seen when qce ~ 1.0. The detailed source types in region A, B and 
C are listed in Tables 2-4, and summarized below. 

(1) The high-luminosity {Lx > W^^ ergs~^) sources in the 10 < i? < 300 pc region in all models are mainly RLOF 
BH-XRBs with HeMS companions. They all have short orbital period (~ 1 hr) and low-mass (~ 1-3 Mq) companions. The 
system velocity is ~ 10-20 kms^^ for models M1-M4, and slightly higher, ~ 30-100 kms~^ for models M5-M8. 

(2) The low- luminosity {Lx < W^^ ergs~^) sources in the 10 < i? < 300 pc regions in all models are mainly (> 80%) 



NS-XRBs with HeMS companions (from 
wind-capture. The orbital periods are ~ 1 



50% in model M4 to ~ 85% in model M7), transferring mass through RLOF or 
- 20 hr for models Ml and M3-M8, and can reach ~ 10^ hr for model M2. The 
1 - 4M(7) for models Ml and M4, ~ 1 _ 20Mpt for models M5-M8 and can reach ~ 60Mo for models 



companion masses are 
M2 and M3. 

(3) The XRBs in the 300 < R < 1000 pc regions in all models all have relatively low- mass (~ 1-4 Mq) companions. Their 
orbital periods are about several hours, and their system velocities are ~ 150-300 kms~^. For models M1-M4, only NS-XRBs 
can be found in this region, with both MS and HeMS donors. For models M5-M8, there are nearly comparable number of 
BH-XRBs and NS-XRBs, with mainly MS donors. The BH-XRBs in these models are from AIC of NS systems. 

Our results indicate that the IMFs of the primary and the secondary stars, and the kick velocity can affect the formation 
and evolution of massive XRBs, but not so significant to constrain models by comparing with the observed Lx vs. R relation. 
However, changing the CE parameter acE can make significant differences in stellar components that construct the correlation. 
We suggest the remarkable different types of XRB in the 300 < R < 1000 pc region may provide an interesting clue to constrain 

the CE parameter. 

Common envelop e evolution has long been the subject s of studie s , from early semi-empiri c al est imate s by Ostrikei (Il975 ) 
and Paczviiskil (j 19761 1 . followed by series of works bv e.g.. IWebbinkI lil992l) . ISandauist etal] (|200d l. and iTaam fc SandauistI 



1 2OO0I 1 ■ to recent 3-D hydrodynamical simulations by Richer fc'^aml (j2008h . However, reliable value of qce is still uncertain 



due to l ack of understanding of the processes involved. P revious semi- empirical estimates suggested that qqe is in the range 
0.6-1.0 (jlben fc Tutukovll 19891 : Frutukov fc Yungelonlll993l ) and 0.3-0.6 jLivio fc Sokeijl 19881 : Ffaam fc Bodenheimei]|l989l ). but 
the value was later suggested to be greater than unity if additional energy source is considered, such as thermal and ionization 
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energy of the envelope (IHan et al. 1995) or possib ly nuclear energy generated in the region close to the in-spirahng companion 
l Taamlll99'j 'l and others fsee llben fc Livioll 19931 . :'or more). In som e evolutionary calculations, high value of acE (~ 3) was 
required in order to model the formation of low- mass BH binaries (IPortegiGS Zwart et al."l997'), low- and intermediate-mass 
BH-XRBs (Kalogcra 1999) and low-mass, short-period BH binaries (Podsiadlowski ct al. 2003). Note in our calculations, much 
more BH-XRBs can be produced in the acE ~ 3.0 cases than in the acE = 1.0 cases, which is consistent with their findings. 

Our results are subject to some uncertainties and simplified treatments. For example, in our calculations, we only consider 
the primordial binaries. However, dynamical interactions in star clusters can produce new XRBs and modify the IMFs of the 
primordial binary stars, though the proportion may be small, as estimated above. Another simplified treatment is taken for 
the CE evolution. Besides acE, the A-parameter, which describes the binding energy between the envelope and the core of 
the donor star, can also affect the CE evolution. It depends on the structure of the stellar envelope, and consequently on 
the evolutionary stage of the star. In our calculation, we adopt a constant (i.e., A — 0.5) disregarding the evolution of the 
donor star at the onset of the mass transfer process, though the value of A changes remarkably. For example, typical value 
of A was found to be between 0.2 an d 0.8, but can b e greater than 5 for the asymptotic giant branch of lower-mass stars 
((Pewi fc TaurislboOOl ). Recent work bv lXu fc Lil l|2O10l ) showed that for more massive stars (i.e , > lOM©), the value of A may 
be low as 0.1. So incorporating A as a function of stellar radius will help model the Lx vs. R relation more precisely, and put 
more realistic constraints on the model parameters, though it is beyond the scope of this paper. 



5 ACKNOWLEDGMENTS 

We are grateful to the referee for his helpful comments and suggestions that greatly improved the paper. We also thank Xi-wei 
Liu and Xiao-jie Xu for useful discussions. This work was supported by the Natural Science Foundation of China (under 
grant number 10873008) and the National Basic Research Program of China (973 Program 2009CB824800). ZYZ was also 
supported by the Jiangsu Project Innovation for PhD candidates (0201001504). 



REFERENCES 

Abt H. A., 1983, ARA&A, 21, 343 

Bailyn C. D., Grindlay J. E., 1987, ApJ, 316, L25 

Ballero S., Matteucci F, Origha L., Rich R., 2007, A&A, 467, 123 

Bassa C. G., van Kerkwijk M. H., Koester D., Verbunt F., 2006, A&A, 456, 295 

Belczynski K., Kalogera V., Bulik T., 2002, ApJ, 572, 407 

Belczynski K., Taam R. E., 2003, ApJ, 616, 1159 

Chen W., Shrader C. R., Livio M., 1997, ApJ, 491, 312 

Dalton W. W., Sarazin C. L., 1995, ApJ, 448, 369 

Davies M. B., Hansen B. M. S., 1998, MNRAS, 301, 15 

Dewi J.D.M., Tauris T.M., 2000, A&A, 360, 1043 

Fabbiano G., 1989, ARA&A, 27, 87 

Fryer C. L., Woosley S. E., Hartmann D. H., 1999, ApJ, 526, 152 

Fryer C. L., 1999, ApJ, 522, 413 

Fryer C, Kalogera V., 2001, ApJ, 554, 548 

Ho L. C, Filippenko A. V., 1996a, ApJ, 466, L83 

Ho L. C, Filippenko A. V., 1996b, ApJ, 472, 600 

Han Z., Podsiadlowski P., Eggleton P. P., 1995, MNRAS, 270, 121 

Hansen B., Phinney E., 1997, MNRAS, 291, 569 

Hobbs C, Lorimer D. R., Lyne A. C, Kramer M., 2005, MNRAS, 360, 963 
Hogeveen S. J., 1990, Ap&SS, 173, 315 

Hurley J. R., Pols O. R., Tout C. A., 2000, MNRAS, 315, 543 
Hurley J. R., Tout C. A., Pols O. R., 2002, MNRAS, 329, 897 
Hut P., Verbunt F., 1983, Nature, 301, 587 
Iben I. Jr., Tutukov A. V., 1989, ApJ, 342, 430 
Iben I. Jr, Livio M., 1993, PASP, 105, 1373 

Garcia M. R., Miller J. M., McClintock J. E., King A. R., Orosz J., 2003, ApJ, 591, 388 
Ivanova N., Kalogera V., 2006, ApJ, 636, 985 

Kaaret P., Corbel S., Prestwich A. H., Zezas A., 2003, Science, 299, 365 

Kaaret P., Alonso-Herrero A., GaUagher J. S., Fabbiano G., Zezas A., Rieke M. J., 2004, MNRAS, 348, L28 
Kalogera V., Webbink R. F., 1998, ApJ, 493, 351 



© 2002 RAS, MNRAS 000, ??-?? 



Displacement of XRBs from star clusters 9 



Kalogera V., 1999, ApJ, 521, 723 

Kaper L., van der Meer A., Massive Stars in Interactive Binaries, ASP Conference Series 367. Edited by Nicole St.-Louis 

and Anthony F.J. Moffat., San Francisco: Astronomical Society of the Pacific, 2007., p. 447 
Kiel P. D., Hurley J. R., 2006, MNRAS, 369, 1152 
Kiel P. D., Hurley J. R., 2009, MNRAS, 395, 2326 

King A. R., Davies M. B., Ward M. J., Fabbiano G., Elvis M., 2001, ApJ, 552, L109 
Kobulnicky H. A., Fryer C. L., 2007, ApJ, 670, 747 
Kording E., Falcke H., Markoff S., 2002, A&A, 382, L13 
Kroupa P., Tout C. A., Gilmore G., 1993, MNRAS, 262, 545 
Kulkarni S. R., Hut P., McMillian S., 1993, Nature, 364, 421 

Larson R. B., 2001, in lAU Symp. 200, The Formation of Binary Stars, ed. H. Zinnecker & R. D. Mathieu (San Francisco: 
ASP), 93 

Liu X. W., Li X. D., 2007, ChJA&A, 7, 389 

Livio M., Sokcr N., 1988, ApJ, 329, 764 

Lucy L. B., 2006, A&A, 457, 629 

Lyne A. G., Lorimer D. R., 1994, Nature, 369, 124 

Masscy P. M., Hunter D. A., 1998, ApJ, 493, 180 

McCrady N., Gilbert A. M., Graham J. R., 2003, ApJ, 596, 240 

Nelemans G., Tauris T. M., van den Heuvel E. P. J., 1999, A&A, 352, L87 

Ostrikcr J. P., 1975, paper presented at lAU Symp., 73, The Structure and Evolution of Close Binary Systems 

Paczyhski B., 1976, In: Eggleton P., Mitton S., Whelan J. (eds.) Structure and Evolution in Close Binary Systems. Proc. 

lAU Symp. 73, Reidel, Dordrecht, p. 75 
Paczynski B., 1990, ApJ, 348, 485 
Phinney E. S., Sigurdsson S., 1991, Nature, 349, 220 
Podsiadlowski P., Rappaport S. A., Han Z., 2003, MNRAS, 341, 385 
Portegies Zwart, S. F., Vcrbunt, F., Ergma, E., 1997, A&A, 321, 207 
Portegies Zwart S. F., Makino J., McMillian S. L. W., Hut P., 1999, A&A, 348, 117 
Pourbaix D. et al., 2004, A&A, 424, 727 

Rappaport S. A., Podsiadlowski P., Pfahl E., 2004, MNRAS, 361, 971 

Ricker P. M., Taam R. E., 2008, ApJ, 672, L41 

Sandquist E. L., Taam R. E., Burkert A., 2000, ApJ, 533, 984 

Sigurdsson S., Hernquist L., 1993, Nature, 364, 423 

Smith L. J., Gallagher J. S., 2001, MNRAS, 326, 1027 

Sternberg A., 1998, ApJ, 506, 721 

Taam R., Bodenheimer P., 1989, ApJ, 337, 849 

Taam R. E., 1994, in ASP Conf. Ser., 56, Interacting Binary Stars, ed. A. W. Shafter (San Francisco: ASP), 208 
Taam R., Sandquist E., 2000, ARA&A, 38,113 

Tauris T.M., van den Heuvel E. P. J., 2006, in Compact Stellar X-ray Sources, ed. W. H. G. Lewin & M. van der Klis, 

Cambridge Univ. Press, (Cambridge), 623 
Tutukov A. v., Yungelon L. R., 1993, MNRAS, 260, 675 
van den Heuvel E.P.J. , 1975, ApJ, 198, L109 

van den Heuvel E. P. J., Portegies Zwart S. F., Bhattacharya D., Kaper L., 2000, A&A, 364, 563 
Van Paradijs J., White N., 1995, ApJ, 447, L33 

Van Paradijs J., 1996, ApJ, 464, L139 
Verbunt P., Hut P., 1987, lAUS, 125, 187 

Verbunt F., van den Heuvel E., 1994, X-ray Binaries, Cambridge University Press, p. 457 
Webbink R. F., Rappaport S., Savonije G. J., 1983, ApJ, 270, 678 
Webbink R.F., 1984, ApJ, 277, 355 

Webbink R. F., 1992, in X-Ray Binaries and Recycled Pulsars, van den Heuvel E. P. J., Rappaport S. A., (NATO ASI Ser. 

C, 377; Dordrecht: Kluwer), 269 
White N.E., Van Paradijs J., 1996, ApJ, 473, L25 
Xu X.J., Li X. D., 2010, in preparation 
Zuo Z. Y., Li X. D., Liu X. W., 2008, MNRAS, 387, 121 



© 2002 RAS, MNRAS 000, ??-?? 



10 Zuo & Li 



Table 1. Parameters adopted for each model. Here acE is the CE parameter, q the initial mass ratio, c-jjick the dispersion of kick speed, 
IMF is the initial mass function. 



Model 




P(q) 


fkick 
km/s 


IMF 


Ml 


1.0 




190 


Kroupa 


M2 


1.0 




190 


Ballero 


M3 


1.0 


oc 


190 


Kroupa 


M4 


1.0 




270 


Kroupa 


M5 


3.0 


ocgO 


190 


Kroupa 


M6 


3.0 


ocgO 


190 


Ballero 


M7 


3.0 


(X q^ 


190 


Kroupa 


M8 


3.0 


(xq° 


270 


Kroupa 



Table 2. The detailed types of sources in region A (Lx > 10^* ergs"!, 10 < i? < 300 pc). Here "BH(NS)RLO", "BH(NS)MS" and 
"BH(NS)HeMS" represent RLOF BH(NS)-XRBs, BH(NS)-XRBs with MS companions and BH(NS)-XRBs with HeMS companions, 

respectively. BH% represents the percentage of BH-XRBs in region A, j!^(^]^QaUgrg^gj represents the percentage of high-luminosity (Lx > 

10'^* ergs~^) sources in 10 < R < 300 pc region. 



Model 


BlI'X. 


DIIRLO 


XSRLO 




Dims 


DIIIIrMS 


NSMS 


NSIIc-MS 


BII 






DIl 


BII 


XS 


NS 


Ml 


99 


99 


100 


22 


1 


97 





100 


M2 


85 


99 


100 


49 


11 


88 





100 


M3 


82 


99 


100 


45 





99 





100 


M4 


99 


96 


100 


24 


11 


84 





100 


M5 


97 


99 


100 


64 


9 


88 


95 


5 


M6 


97 


99 


100 


67 


6 


90 


91 


9 


M7 


97 


99 


100 


68 


1 


95 


43 


57 


M8 


98 


99 


100 


60 


5 


92 


96 


4 



Table 3. Same as in Tabic 2 but for sources in region B (lO'^*^ < Lx < 10^* ergs"!, 10 < R < 300 pc). Here ^ 

represents the percontagc of low-luriiiiiositA' (10"*'' < Lx < 10'""^ crgs^^) sonrcos in 10 < i? < 300 pc region. 



Model 


BH% 


BHRLO 


NSRLO 


N(10^'5<Lx<10^*erg/s) 


BHMS 


BHHeMS 


NSMS 


NSHeMS 


BH 


NS 


N(>10!*«erg/s) 


BH 


BH 


NS 


NS 


Ml 


3 


2 


28 


78 


70 


28 


11 


80 


M2 


11 





58 


51 


80 


18 


23 


70 


M3 


8 


34 


58 


55 


52 


41 


8 


76 


M4 


7 


19 


51 


76 


62 


20 


21 


47 


M5 


13 


63 


57 


36 


37 


38 


33 


64 


M6 


20 


46 


57 


33 


29 


48 


27 


70 


M7 


11 


29 


44 


32 


13 


64 


11 


84 


M8 


9 


34 


57 


40 


38 


56 


35 


61 
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Table 4. Same as in Tabic 2 but for sources in region C (Lx > 10^^ ergs ^ , 300 < -R < 1000 pc) . represents the percentage 

of high-luminosity (Lx > 10"^* ergs^^) sources in region C. 



iVLoaei 


r>xl/o 


BHRLO 
BH 


NSRLO 
NS 


N(>10^*erg/s) 
N(>loaBerg/s) 


BHMS 
BH 


BHHeMS 
BH 


NSMS 
NS 


NSHeMS 
NS 


Ml 








82 











37 


44 


M2 








23 











15 


69 


M3 








33 











26 


0" 


M4 








60 











21 


33 


M5 


38 


100 


95 


14 


95 


5 


90 


4 


M6 


59 


100 


99 


19 


96 


4 


95 


1 


M7 


35 


100 


99 


21 


71 


29 


98 





M8 


40 


100 


96 


9 


99 


1 


88 


2 



" naJjed Helium star Hertzsprung Gap: 67% 
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36 I , , , I , , , I , , , I , , , 

1.0 1.5 2.0 2.5 1.0 1.5 2.0 2.5 1.0 1.5 2.0 2.5 1.0 1.5 2.0 2.5 3.0 

Log(R/pc) 

Figure 3. The Lx - R distribution for ALL-XRBs (top), NS-XRBs (middle) and BH-XRBs (bottom), respectively. Prom left to right 
are models M1-M4, respectively. 




36 I . . . Ib . . . I . . . I . . . I 
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Log(R/pc) 

Figure 4. Same as Fig. 3 but for sources in models M5-M8 from left to right, respectively. 
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Figure 5. The Lx — M2, Lx — Porht ^^nd Porb ~ ^2 distributions in the 10 < iJ < 300 pc region from left to right, respectively. From 
top to bottom are models M1-M4, respectively. 
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-2-1012-10123-2-1012 
Log(M2) Log(P„,,/hr) loq{U^) 

Figure 6. Same as Fig. 5 but for sources in models M5-M8 from top to bottom, respectively. 
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Figure 7. The Lx — ^2, Lx — forbi ^orb ~ ^2 distributions in tiie 300 < R < 1000 pc region from left to rigiit, respectively. From top 
to bottom axe models M5-M8, respectively. 



© 2002 RAS, MNRAS 000, ??-?? 



Displacement of XRBs from star clusters 17 




i 



M5 



H — I — I — I — ■ — I — h 



M6 



H — I — I — I — 1 — I — I- 



M7 



H — I — I — I — ■ — I — h 



39 

38 

37 

36 

39 
38 
37 



o 



39 

38 

37 
36 

39 

38 

37 
36 



H — I — I — I — ■ — I — I- 



H — I — I — I — 1 — I — I- 



H — I — I — I — ■ — I — I- 



■1 1 
Log(M2) 



1 2 
Log(Porb/hr) 



o 



H — I — I — I — ■ — I — h 



H — I — I — I — 1 — I — I- 



H — I — I — I — ■ — I — h 



1 1 
LogCMj 



Figure 8. Same as Fig. 7 but for sources in models M5-M8 from top to bottom, respectively. 
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Figure 9. The delay time distributions between SN and the beginning of RLOF for sources in regions A (dash-dot-dotted Une), B (solid 
line) and C (dotted line), respectively. From left to right are models M1-M4 (upper panel) and M5-M8 (bottom panel), respectively. 
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Figure 10. The Porb.SNe ~ -'^2,SNe distributions in regions A, B and C, respectively. Prom left to right are models M1-M4, respectively. 
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Figure 11. Same as Fig. 10 but for sources in models M5-M8 from left to right, respectively. 
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Figure 12. The evolution of Mi, M2, Porbi ^■nd Lx for an example XRB in region C. 
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Figure 13. Same as Fig. 3 but without AIC of accreting NSs. 




Figure 14. Same as Fig. 4 but without AIC of axjcreting NSs.. 
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Figure 15. The evolution of tlie Lx — R distribution for models Ml (top panel) and M5 (bottom panel), respectively. From left to right 
are for the ages of 20 Myr (left), 30 Myr (middle) and 40 Myr (right), respectively. 
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